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ABSTRACT 

In this paper, we assess the constraints on the evolutionary models of young low-mass objects that 
are provided by the measurements of the companion AB Dor C by Close and coworkers and by a 
new comparison of model-derived initial mass functions (IMFs) of star-forming regions to the well- 
calibrated IMF of the solar neighborhood. After performing an independent analysis of all of the 
imaging and spectroscopic data for AB Dor C that were obtained by Close, we find that AB Dor C 
(which has no methane) is not detected at a significant level (signal-to-noise^ 1.2) in the simultaneous 
differential images (SDI) when one narrow-band image is subtracted from another, but that it does 
appear in the individual SDI frames as well as the images at J, 77, and Kg- Our broad band photometry 
for AB Dor C is consistent with that of Close. However, the photometric uncertainties that we measure 
are larger than those derived by Close; our uncertainties are consistent with those measured in other 
studies using the same adaptive optics system. Using the age of r = 75-150 Myr recently estimated for 
AB Dor by Luhman, Stauffer, and Mamajek, the luminosity predicted by the models of Chabrier and 
Baraffe is consistent with the value that we estimate from the photometry for AB Dor C. We measure 
a spectral type of M6±l from the K-hand spectrum of AB Dor C, which is earlier than the value of 
M8±l reported by Close and is consistent with the model predictions when a dwarf temperature scale 
is adopted. In a test of these evolutionary models at much younger ages, we show that the low-mass 
IMFs that they produce for star-forming regions are similar to the IMF of the solar neighborhood. 
If the masses of the low-mass stars and brown dwarfs in these IMFs of star-forming regions were 
underestimated by a factor of two as suggested by Close, then the IMF characterizing the current 
generation of Galactic star formation would have to be radically different from the IMF of the solar 
neighborhood. 

Subject headings: instrumentation: adaptive optics — binaries: visual — stars: formation — stars: 
low-mass, brown dwarfs — stars: pre-main sequence 



1. INTRODUCTION 

For decades, theoretical evolutionary models have been 
an essential tool for interpreting the positions of stel- 
lar populations on the Hertzsprung-Russell diagram. In 
particular, studies of star and planet formation have re- 
lied on these models to provide estimates of masses and 
ages of members of star- forming regions (r < 10 Myr). 
Because the uncertainties in the models are largest at 
such early stages (Baraffe ct al. 2002), observational 
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masses below 0.3 Mp^, fewer tests of the mass esti- 
mate s are availablejlM^h^n^e^^j 120041: iReiners et al.l 
EoHS)- Recently, iClose et all l|2005r partially addressed 
this problem by measuring a dynamical mass for a low- 

^ Based on observations made with ESO Telescopes at the 
Paranal Observatories under program ID 60.A-9026. This pub- 
lication makes use of data products from the Two Micron All Sky 
Survey, which is a joint project of the University of Massachusetts 
and the Infrared Processing and Analysis Center /California Insti- 
tute of Technology, funded by the National Aeronautics and Space 
Administration and the National Science Foundation. 

^ Department of Astronomy and Astrophysics, The Pennsyl- 
vania State University, University Park, PA 16802, USA; kluh- 
man@astro.psu.edu. 

^ Steward Observatory, The University of Arizona, Tucson, AZ 
85721, USA; potter@as.arizona.edu. 



mass companion, AB Dor C, that is much older than 
star-forming regions but still above the main sequence. 
This object was first detected indirectly bv lGuirado et alJ 
(1997) through astrometry of AB Dor, which is a 
well-studied nearby young K-type star {d = 14.9 pc). 
[Close et al. (2005) resolved AB Dor C from the primary 
with adaptive optics (AO) images and u sed a combina- 
tion of these data and the astrometry from lGuirado et all 
(1997) to arrive at a mass of M = 0.09 ± 0.005 Mq. 
After adopting an age of 50 Myr for the AB Dor sys- 
te m, they found that the theoretical evolutionary models 
of Ch abrier et al.1 l)2000D overestimated the near-infrared 
(IR) flu xes of AB Dor C b y roughly one magnitude. As 
a result. IClose et all l)2005D concluded that these models 
underestimate the masses of young low-mass objects by 
a factor of two and that many of the objects previously 
identified as brown dwarfs in star-forming regions and 

open clusters are instead low-mass st ars. 

The observations of AB Dor C bv IClose et all l|2005fl 
(hereafter C05) potentially provide an important new 
test of the theoretical evolutionary models of low-mass 
objects. However, the results of this test hinge on reliable 
measurements of t he mass, age, luminosity, and tempera- 
ture of AB Dor C. lLuhman. Stauffer. k, MamaieS l|2005fl 
recently reexamined the age constraints of the AB Dor 
system, concluding that it has the same age as the 
Pleiades open cluster (r = 100-125 Myr), which removed 
the discrepancy between the luminosity r eported by C05 
and the value predicted by the models of iChabrier et all 
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l)2000tl . In this paper, we reexamine other steUar prop- 
erties of AB Dor C, namely its luminosity and tem- 
perature, and describe the constraints these data pro- 
vide for evolutionary models (§ EJ- In addition, we 
present a second test of the models at low masses, but 
now at much younger ages (r ~ 1 Myr), in which we 
compare the low-mass initial mass function (IMF) for 
nearby star-for ming re gions as derived with the models 
of Baraffc ct aQ l|1998|) and^^habrier et al. ( 2000) to the 
relatively accurate and well-calibrated IMF of the solar 
neighborhood under the assumption that the true IMFs 
of these populations are the same (§|S1)- 

2. AB DOR C 

2.1. Photometry 

Given the large flux ratio and small angular separation 
of AB Dor A and C, a robust photometric measurement 
of the latter with AO images requires careful considera- 
tion of the possible errors inherent to data of this kind. 
For these reasons, we have analyzed the data published 
in COS, which are now available to the public at the data 
archive of the European Southern Observatory (ESO). 

2.1.1. SDI Images 

The data with which COS reported the first direct de- 
tection of AB Dor C were obtained during the commis- 
sioning of the simu ltaneous differential imaging mode 
(SDI, Lcn zen et aLl '2004: Billcr ct al. 2004) of the High- 
resolution Near-IR Camera (CONICA, Lenz en et al.l 
[1998), which was used in conjunction with the Nasmyth 
Adaptive Optics System (NAOS, Rousset et al. 2000) on 
ESO's Very Large Telescope (VLT). T he SDI technique, 
first ilhi strated bvlRacine et alJ lll999D and further devel- 
oped bvlMar ois et all l)2000ll2003ll200S D. uses the relative 
spectral contrast between a primary and a companion to 
suppress the contrast-limiting speckle noise. It has been 
shown that significant improvements in speckle suppres- 
sion are gained through the SDI technique. However, 
non-common path errors, differential chromatic aberra- 
tions, alignment-sensitive ghost images, flat-field miscal- 
ibration, and detector read noise can significantly exceed 
the photon noise. Although these data are not directly 
relevant to testing the theoretical evolutionary models, 
we examine them in this section because of the impor- 
tance of testing the abilities of new observational tech- 
niques like SDI. 

The SDI observations from COS were obtained on the 
night of 2004 February 1 and consisted of a set of im- 
ages at five dither positions followed by a second set of 
dithered images with the instrument rotator offset by 
33°. At a given position, four SDI images were obtained 
through narrow-band filters, one at 1.S7S /xm, one at 
1.600 /im, and two at 1.62S fim. After removing bad 
pixels and dividing by flat field images, we created im- 
ages that were spatially scaled by the ratio of central fil- 
ter wavelengths to make the diffraction and speckle pat- 
tern for all images coincident in image coordinates. The 
images were optimally aligned to subpixel accuracies us- 
ing the flux conserving interpolation routine interpolate 
within IDL. After the subpixel alignment, the individ- 
ual images were optimally flux scaled to minimize the 
residual counts in the non-saturated regions of the PSF 
halo in the subtracted image. The images were then spa- 



tially flltered as in iMasciadri et"al] l|2005j) . The six dif- 
ferent non-redundant subtraction combinations possible 
between the four images were created using these optimal 
subtraction routines for each of the 10 dither positions. 
Because the plate scale is independent of the wavelength 
in the SDI mode of CONICA, when the images are spa- 
tially de-magnified to match the diffraction pattern of the 
PSFs with different wavelengths, the positions of a com- 
panion in different filters becomes radially misaligned to 
a degree proportional to the demagnification factor and 
the distance from the center of the primary PSF. As a 
result, the signature of a real (non-methane) compan- 
ion using the SDI technique is not just a fainter version 
of the primary PSF but rather a subtraction between 
two radially misaligned PSFs, as illustrated in Figure ^ 
This effect essentially improves the sensitivity to the de- 
tection of companions by producing a more unique and 
detectable signal even when the companion has no intrin- 
sic spectral contrast between the narrow-band filters. 

At the position of AB Dor C reported by COS, we find 
a point source appearing in all four of the unsubtracted 
images of the SDI data at a signal-to- noise of ~ 3. How- 
ever, this feature overlaps with the diffraction spider arm, 
making the differential imaging crucial in confirming it 
is as a real companion with these data alone. Indeed, 
after the optimal alignment and subtraction routine was 
performed through the dataset, we found no obvious sig- 
nature of a companion at this location in any iteration 
of the optimally subtracted images. Although not done 
by COS, we performed an additional step to remove pos- 
sible fixed noise from the reduced differential images by 
subtracting the reduced data sets taken at the differ- 
ent instrument rotator angles. This would produce a 
negative version of the companion's signal at the same 
radius as the companion, but separated in position an- 
gle equal to the instrument rotator angle change (33°). 
From this analysis, we again found no obvious signa- 
ture of a co mpanion. A similar t ype of reduction was 
performed bv iNielsen et alJ l|2005D . resulting in the im- 
ages in the first panel of their Figure 4. In contrast to 
their claims, we see no significant detection of a com- 
panion in those data. To quantify the detection limits, 
a noise analysis was performed using the five dithered 
images for each rotation angle for each of the six sub- 
traction combinations. Because the central five pixels 
surrounding AB Dor A were saturated, a PSF for use 
with artificial companions could not be measured from 
these data. However, prior to that sequence of images, 
COS obtained a single short exposure in which AB Dor A 
was not saturated. We adopt this PSF of AB Dor A 
for artificial companions in the following analysis. This 
use of a PSF from a short exposure tends to overesti- 
mate the Strehl ratio of the longer exposure, which re- 
sults in a detection limit that is overly optimistic. Based 
on the spectra of late-M objects at a variety of ages (e.g., 
Cushing et al. 200S), the maximum ffux ratio expected 
for AB Dor C between the 1.62S and 1.S7S /im narrow- 
band filters is i^(1.62S)/F(l.S7S) ~ 1.1. Therefore, we 
simulated AB Dor C with an artificial companion that 
was 125 and 112 times fainter than AB Dor A for the 
images at 1.625 and 1.S7S /im, respectively. 

By combining the signal from the artificial companion 
with the standard deviation of the total flux in bins of 
linear size similar to the diffraction limit (^ 4 pixels) 
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around a ring at the radius of AB Dor C reported by 
C05 (~ 9 pixels), we arrived at an expected signal-to- 
noise ratio of 1.2 for AB Dor C. We checked the validity 
of this measurement by inserting the artificial compan- 
ion into the raw data at the separation of 0'.'16 reported 
by COS, repeating the reduction procedures, and visu- 
ally inspecting the resulting images. For comparison, we 
also included companions at larger separations of 0'.'4, 
0'.'64, and 0'.'88. The differenced SDI images containing 
these companions are shown in Figure ^ Although a 
companion of this contrast would easily be detected at 
separations beyond O'.'S, it is not reliably detected at the 
radius of (/.'IG reported by COS for AB Dor C. Thus, we 
cannot reproduce the putative detection of a companion 
shown in Figure lb from COS. Their image also does not 
exhibit the expected plus/minus asymmetry expected for 
SDI data of this kind. Nor is that source evident in the 
difference between two instrum ent rotator angle s in the 
first panel of Figure 4 fr om Ni elsen et alJ l)200SD . How- 
ever, NielsciTe^^ l)200S() does show a significant detec- 
tion of AB Dor C in newer SDI data obtained after COS 
when the companion had moved to a larger separation, 
which is consistent with the detection limits that we have 
measured. 

2.1.2. Broad-band Images 

In addition to the SDI data, COS obtained near-IR 
broad-band images of AB Dor with CONICA and NAOS 
(NACO). For AB Dor C, COS reported two a errors of 
+0.19 and -0.24 at J, +0.13 and -O.IS at H, +0.12 and 
-O.IS at Kg. These errors correspond to signal-to-noise 
ratios of 10, 16, and 11 in the three bands, respectively, 
which are substantially better than achieved in previ- 
ous companion search es with NACO l|Biller et al.ll2003 
iMasciadri et al.ll2005j) . To investigate this apparent dis- 
crepancy, in this section we present a detailed, explicit 
analysis of the photometric errors for AB Dor C in the 
data obtained by COS, including the possible systematic 
errors inherent to high contrast photometric measure- 
ments with AO observations. We also consider the un- 
published images in a narrow-band filter at 3.74 fim that 
were obtained of AB Dor with NACO. 

On the night of 2004 February 4, COS obtained IS im- 
ages with exposure times of 10 s among S dither positions 
for each of the J, H, and Kg filters and three images with 
exposure times of 3.4S s (10 coadds of 0.345 sec). A sim- 
ilar set of exposures was obtained through the NB3.74 
filter, but with exposure times of 3.6 s (10 coadds of 
0.36 sec). After flat fielding, dark subtraction, and re- 
moval of bad pixels, we constructed a sky frame in a 
given filter from a median combination of unaligned im- 
ages and subtracted it from each image in that filter. We 
aligned the resulting images to subpixel accuracy using 
the same IDL interpolation routine as in the SDI data 
reduction. In the long exposures, pixels within a radius 
of ~ 0'.'07 from AB Dor A were saturated. These images 
were aligned using the unsaturated wings of AB Dor A. 
We also generated a set of images that were spatially 
filtered by subtracting from each image the same image 
convolved with a Gaussian with a width roughly equal to 
the FWHM for each wavelength (Masciadri ct al. 200S). 
In the resulting images, the low frequency PSF variations 
are suppressed. The final nonfiltered and filtered images 
for J, H, and Kg are shown in Figure |21 



The companion AB Dor C is apparent in the individual 
reduced 10 s images at J, H, and Kg, but not NB3.74. A 
detection at NB3.74 is not expected given the total expo- 
sure time in that filter. Although other fixed speckles and 
diffraction artifacts are present in the halo of AB Dor A 
at comparable brightness, AB Dor C is the only "speckle" 
that remains at a constant position in the J, H, and Kg 
images. If it were a fixed diffraction speckle, the sepa- 
ration would be proportional to the central wavelength 
of the filter. The fact that it is in the same position in 
the three broad-band filters as well as the unsubtracted 
SDI images and is in all the different dither positions 
demonstrates the reality of the companion. 

To measure photometry for AB Dor C from the J, H, 
and Kg images, we used the shorter, unsaturated expo- 
sures of AB Dor A to construct a PSF for each filter 
and calibrated it with photometry of AB Dor A from 
2MASS. We confirmed the accuracy of this calibration 
by extracting photometry for AB Dor Ba-|-Bb, which ap- 
peared in some of the long exposures, and comparing it 
to the 2MASS measurements. In each filter, the cali- 
brated PSF was scaled and subtracted from the image 
at the position of AB Dor C with a scale factor opti- 
mized such that the residual at the location of AB Dor C 
matched the light level of the adjoining regions. The er- 
rors in these photometric measurements were computed 
from the standard deviation of the brightness at the po- 
sition of AB Dor C in an area with a width equal to 
the PSF's FWHM through the IS images. In addition, 
we calculated the standard deviation of flux in FWHM- 
sized bins at all position angles at a constant radius from 
AB Dor A as done in the con trast sensitivity estima- 
tions of IMasciadri et all l)200SD . This analysis resulted 
in errors slightly larger than those calculated from the 
variance through a data set at the position of AB Dor C; 
specifically, the e rrors are 12%, 17%, and 7% larger using 
the method from lMasciadri et all pOOSi) for the J, H and 
Kg bands, respectively. This is likely due to the fact that 
the diffraction spikes from the secondary support create 
intrinsic variability around the PSF halo at a constant 
radius. For this reason, we use the lower error values 
calculated as the standard deviation of the flux through 
the data set at the position of AB Dor C. The photom- 
etry and errors produced by this analysis are presented 
in Tabled We also list the measurements from COS for 
comparison. The photometric errors for AB Dor C that 
we derive from COS's images are significantly larger than 
those reported by COS. In compar i son, t he nonfiltered 
NACO data from IMasciadri et all l)200Sj) exhibited an 
average sensitivity of S a at Am ~ 4.S ± O.S and 0'.'16, 
which corresponds to '--^ 3 cr for an object at the contrast 
of AB Dor C and thus is consistent with the photometric 
errors we measure from the no nfiltered images. Note t hat 
although the sensitivities from lMasciadri et al.l l)200S() for 
filtered images are better than those from the nonfiltered 
data, they represent only the significance of a detection 
rather than bona fide photometric errors because they 
do not include the systematic errors introduced by the 
filtering process. 

In addition to the photometric errors listed in Table [H 
the noise associated with the variability of the reference 
PSF must also be taken into consideration. The refer- 
ence PSFs are from the short exposures of AB Dor A 
(0.34S s X 10 coadds x 3), which were measured at a 
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different time than the longer exposures that detected 
AB Dor C (10 s x 1 coadd x 15). Because the total ex- 
posure time among the short exposures in a given band 
was comparable to just one long exposure, we cannot di- 
rectly measure the temporal variability of the PSFs used 
in the photometric analysis. However, the NB3.74 data, 
which were obtained on the same night as the long expo- 
sures, encompassed a longer time baseline. These data 
consisted of a continuous series of 15 frames of 0.36 s 
X 10 coadds. The combination of three consecutive im- 
ages within this series has a similar total clock time as 
one long exposure. Therefore, from the series of 15 im- 
ages we were able to construct 5 independent unsatu- 
rated PSFs that were measured on the same timescale as 
the long exposures. The variability of the NB3.74 PSF 
can then be translated into variabilities at other wave- 
lengths using a relation between the Strehl ratio (S), the 
RMS wavefront error (RMS), and the observing wave- 
length X: S = exp{-{2TTRMS/\f). Defining the frac- 
tional error in the Strehl ratio to be / = [S -\- as)/S, 
the fractional error relationship between two wavelengths 
then can be arrived at through algebraic manipulation: 
/(Ai) = /(A2)'^2/^?). The standard deviation of the 
Strehl ratio of the 5 NB3.74 images was measured to be 
f{^NB3.74:) = 1.028. To match the aperture used in the 
photometric analysis, we estimated the changes in the 
Strehl ratio using the total counts within the FWHM 
core of the PSF. The /(J), /(i?), and f{K,) fractional 
Strehl ratio errors are found to be 1.286, 1.155, and 1.087 
respectively. These errors were propagated with those 
found from the halo noise at the radius of AB Dor C and 
are presented in Table 

2.2. Spectral Type 

C05 used the spectroscopic mode of NACO to obtain 
near-IR spectra of AB Dor C on the night of 2005 Febru- 
ary 4. The NACO system was configured to provide 
a slit width of 0'.'086, a spectral coverage of 2-2.5 /im, 
and a resolution oi R = A/AA ^ 1500. Eight dithered 
one-minute exposures were obtained with the slit aligned 
along AB Dor and the position of AB Dor C measured by 
C05, and another set of eight exposures were taken with 
the slit rotated by 180°. The primary w as saturated in 
these data. C05 and lNielsen et al.l l)2005() described some 
of the details of their processing of these data. 

We have retrieved the NACO spectroscopic data for 
AB Dor C from the ESO archive and have performed an 
independent analysis, which we now describe. After di- 
viding by flat field images, we straightened the images to 
precisely align the dispersion and spatial directions with 
the image coordinate axes. To detect AB Dor C within 
the wing of the PSF of AB Dor, we differenced spectral 
images at opposite instrument rotator angles. To account 
for changes in the PSF, prior to this differencing we mea- 
sured the fluxes as a function of wavelength between radii 
of 8.5 and 12.5 pixels on the left and right sides of AB Dor 
for each of the eight images at each of the two rotator 
angles, which we refer to as F\{i,j,k) where i =left or 
right, j — 1-8, and = or 180 (degrees). We then com- 
puted rx{left) = Fx{left,jl,0)/Fx{left,j2,180) and 
rxiright)^ = Fx{right, jl,0)/ Fx{right,j2, 180) for aU 
combinations of jl and j2. These ratios measured the 
change in the PSF on each side of AB Dor between a 
given pair of exposures. To achieve a reliable PSF sub- 



traction, the PSF change should be the same on each side 
of AB Dor. Therefore, we considered only pair subtrac- 
tions in which rx (left) jrx (right) was approximately con- 
stant with wavelength. For each of these retained pairs, 
we fit a function to rx{left) + rx{right) that was linear 
with wavelength, multiplied the exposure at the second 
rotator angle (fc — 180) by this fit, and subtracted the 
resulting image from the image at the first angle (fc = 0). 
This process produced 9 good detections of AB Dor C 
out of the 16 that were possible. We extracted these 9 
spectra, combined them, and divided by the spectrum of 
a telluric standard star, HD 34286 (G3V). The resulting 
spectrum was multiplied by the solar spectrum to cor- 
rect for the spectral slope and absorption features of this 
standard. Because the 9 individual spectra of AB Dor C 
exhibited significant differences in their spectral slopes, 
the slope of the combined spectrum probably was not 
accurate. We removed the slope of the spectrum by di- 
viding by a polynomial fit to the continuum. The final 
spectrum of AB Dor C is shown in Figure 13 

For the spectral classification of AB Dor C, we use 
the dwarf standards obtaine d with SpeX at the N ASA 
Infrared Telescope Facility bv lCushing et al.l l|2005D . The 
resolution of the SpeX data is i? ~ 2000, which is higher 
than that of the NACO spectrum of AB Dor C (i? ~ 
1500). To enable a reliable comparison of the spectral 
features in these spectra, we smoothed the SpeX data to 
the same resolution as the NACO data. To do this, we 
compared the Ar lamp spectra obtained with the SpeX 
and NACO spectra and identified a Gaussian function 
convolution for the SpeX data that produced the same 
FWHMs as in the NACO spectra. We then applied this 
Gaussian function to the spectra of the M dwarfs from 
CushinE et al. (2005). The resulting SpeX data of Gl 213 
(M4V)," Gl 51 (M5V), Gl 406 (M6V), vB 8 (M7V), and 
vB 10 (M8V) are plotted with the NACO spectrum of 
AB Dor C in Figure |31 The slopes of these standard 
spectra have been removed in the same manner as for 
AB Dor C. 

We now measure the spectral type of AB Dor C by 
comparing the strengths of the Na I doublet, Ca I triplet, 
and CO band heads in its spectrum to those exhibited 
by the dwarf standards in Figure 01 The Na I absorp- 
tion in AB Dor C is weaker than that of M5V and M6V 
and stronger than that of M4V, M7V, and M8V. Be- 
cause this doublet i s weaker in pre-mai n-sequence ob- 
jects than in dwarfs l)Luhman et al. 1998a'l. the strength 
observed for AB Dor C is consistent with types of M5- 
M6, but not M4, M7, or M8. In terms of both Ca I 
and CO, AB Dor C is best matched by M6V. Com- 
pared to M4V and M5V, the lower surface gravity of 
AB Dor C co uld explain its stronger CO, but not the 
weaker Ca I l|Luhman fc Riek^ 119981) . Meanwhile, the 
opposite is true for M7V and M8V. In particular, CO 
is significantly weaker in AB Dor C than in M8V, and 
the inconsistency is even larger when the gravity differ- 
ences are considered. Taken together, the strengths of 
Na I, Ca I, and CO indicate a spectral type of M6±l for 
AB Dor C. 

Based on the spectrum of AB Dor C t hat they mea- 
sured from the NACO data, C05 and iNielsen et alJ 
l|2005D derived a spectral type of M8±l for AB Dor C. 
We comment on a few aspects of their classification. In 
agreement with our results in Figure 01 C05 found that 
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AB Dor C exhibited stronger Na I than vB 10 (M8V). 
However, they attributed this difference to the lower sur- 
face gravity of AB Dor C , even though Na I transit ions in 
gener al (Martf n"et al.ll9 96: Luliman 1999; Gorlova ~et al.l 
|2£)03; McGovcrn et al. 2004) and this X-band doublet 
specifically iL uhman et al. ,1998a l are known to become 
weaker, not stronger, with decreasing surface gravity. 
The stronger Na I absorption in the spectrum from 
COS relative to M8V can only be explained by a spec- 
tral type that is earlier than M8V since the strength of 
this feature reac hes a maximum at M6V and decreases 
with later types ijLuhman fc Riekd ll998t iC ushing et al.l 
UnoBI). In addition to dwarfs. iNielsen et all l|2005D used 
young objects in Upper Scorpius as standards. Com- 
pared to a spectrum of a M7 member of that association 
i|Gorlova et al. 2 0Q3), AB D or C exhibited stronger Na I 
and CO, which INielsen et'a l. (2005) cited as evidence 
for a spectral type later than M7. However, as already 
noted, stronger Na I would imply a type earlier than M7, 
not later. Meanwhile, the difference in CO is likely due 
to the fact that th e spectral resolution o f the data for Up- 
per Scorpius from iGorlova et aP l)2003j) was much lower 
than the resolution for AB Dor C. 

C05 also found good agreement between the broad 
spectral shapes (i.e., steam absorption) of their spectrum 
of AB Dor C and a spectrum of vB 10 from Wilking et al. 
l)1999l) . which seemed to support an M8 classification for 
the former. However, the validity of this comparison is 
questionable for the following reasons. First, C05 cor- 
rected for telluric absorption in their data by dividing by 
the spectrum of a solar- type star. They multiplied this 
result by the solar spectrum to remove the spectral slope 
and absorption featu res intrinsic to the telluric stan- 
dard. In comparison, I Wilking et al.l l)1999r ) also divided 
their spectrum of vB 10 by a telluric standard (an AO 
star), but did not attempt to remove the intrinsic spec- 
trum of that s t andar d. Thus, the procedures of C05 and 
iWilking et al.l l)1999() should produce spectral slopes that 
differ systematically, which would preclude a meaningful 
comparison of the steam bands in the resulting spectra of 
AB Dor C and vB 10. In addition, Nielsen ct al. (2005) 
stated that they (and presumably C05) were unable to 
"preserve the continuum" of AB Dor C during their data 
reduction, which also would obviate the use of steam in 
the spectral classification. 

It is useful to consider the spectral type of AB Dor C in 
the context of the low-mass Pleiades members PPL15 A 
and B, which have well-constrained masses and spectral 
types and have similar ages as AB Dor C. According to 
the spectral type of M8±l reported by C05, AB Dor C 
should be coole r than PPL15 A and B (M6 and M7; 
iBasri fc Martfnl[l99 9). However, the dynamical mass of 
AB Dor C (M = 0.09 ± 0.005 Mq, C05) is equal to 
or greater than the masses of PPL15 A and B, each of 
which have a firm upper mass limit of 0.1 Mp, based 
on the presence of L i (Basri. Marcv. fc Grah am 199^ 
IBasri fc Martini Il999|) '* an d are probably brown dwarfs 
based on their binary data ([Basri fc Martini 1999(1 . Thus, 
these relative spectral types and masses of AB Dor C and 
PPL15 are incompatible. Meanwhile, the spectral type 

* This mass limit is based on theoretical relationships between 
Li abundance, mass, and age, which have been sho wn to be robust 
IBildsten et al 11199^ rEurke. Pinsonneault. fc ijilTstaOO-jl) . 



of M6±l that we measure for AB Dor C is perfectly con- 
sistent with PPL 15. 

2.3. Age 

ILuhman. Stauffer. fc Mamaiek! ("2005") recently derived 
an age for AB Dor from a color-magnitude diagram and 
the kinematics of the AB Dor moving group. In a dia- 
gram of Mk versus V — Kg, they found that the AB Dor 
group and the Pleiades were approximately coeval (r = 
100-125 Myr, IMevnet et aLllT993t IStauffer et aD IT99^. 
This result was supported by the kinematic analy- 
sis, which suggested that the two populations orig- 
inated in th e sam e large scale star- formation event. 
INielsen et all l|20(^5D briefly addressed the age of AB Do r 
and the results of Luhman. Stauffer. & Mamaiek ( 2005^. 
First, they quoted Luhman, St auffer. fc Mamaiek (,20C)5i) 
as finding that the AB Dor group is brighter by 
0.1 mag than the Pleiades in Mk versus V — 
Ks- However, that was not the case. Instead, 
ILuhman. Stauffer. fc Mamaiekl (|2005D reported an offset 
of 0-0.1 mag in a visual comparison of the histograms of 
the Mk offsets between the observed positions of stars 
in each population and a fit to the lower envelope of the 
Pleiades sequenc e. In a more quantitative and defini- 
tive comparison, Luhm an. Stauffer. fc Mamaiekl l(2005f l 
found that the mean offsets of the two populations were 
indistinguishable, in dicating no detectable age difference. 
INielsen et al.l l(2005D then presented a comparison of the 
Pleiades and the AB Dor group sequences in Mk versus 
J — Kg and claimed to find an offset of 0.15 mag in Mk 
between the sequences at J — Kg > 0.4. However, a di- 
agram of that kind is a poor choice for measuring ages 
and comparing sequences because the dynamic range of 
J — Kg is very small, nearly an order of magnitude lower 
than that oiV — Kg. In fact, the sequences of stellar pop- 
ulations are vertical in Mk versus J — Kg for low-mass 
stars later than K7, as illustrated in Figure^ As a result, 
the stars that are normally most valuable for measuring 
ages because of their large displacements above the main 
sequence are made useless by plotting them in Mk ver- 
sus J — Kg. Finally, even if the AB Dor group and the 
Pleiades are compared in Mk versus J — Kg , as we have 
done in Figure 0] their sequences are still very similar. 
Following the procedure used for Mk versus V — Kghy 
ILuhman. Stauffer. fc Mamaiekl ((20051) , we compared the 
mean offsets of the stars in the AB Dor group and the 
Pleiades from the lower envelope of the Pleiades. The dif- 
ference between these mean offsets ranges between -0.05 
and 0.05 depending on the exact range of colors that 
is considered. Thus, the mean offsets for the AB Dor 
group and the Pleiades are indistinguishable in Mk ver- 
sus J —Kg, in agreement with the results fr om Mk versus 
V — K g , and we find no basis in the work of INielsen et alJ 
i2005l) for modifying the age estimat e of AB Dor from 
ILuhman. Stauffer. fc Mamaiekl (|2005D . 

2.4. Comparison to Model Predictions 

We now estimate the bolometric luminosity and effec- 
tive temperature of AB Dor C from our broad-band pho- 
tometry and spectral classification and compare them 
to the values predict ed by the evolutionary models of 
IChabrier et"aITl|2000(l . 

We have computed the luminosity by combining our 
corrected H magnitude in Table ^ the distance of 
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AB Dor ijPerrvman et alJll997j) . the g-ba nd bolometric 
corre c tion of M6 f i eld dw arfs (BCh — 2.6. lTinnev et al.l 
119931 iDahn et~al\ l2002j) . and an absolute bolometric 
magnitude of 4.75 for the Sun, arriving at Lboi — 
0.00141° ™°^ Lq. Because the J - H and H ~ Ks col- 
ors of AB Dor C are consistent with those of M6 field 
dwarfs, a luminosity derived from J or Kg produces a 
similar result. In Figure El we compare our measure- 
ments of Mj, Mh, Mk, and Lboi for AB Dor C to the 
values pr edicted by the theoretical evolutionary mod- 
els of iChabrier et alJ ll200Q["). For AB Dor C , wc adopt 
an age of 75-150 Mvr ijLuhman. Stauffer. fc Mamaie'3 
I2005|l . Using this age, the luminosity measured by C05 
agrees well with the model p r edicti ons, as shown by 
iLuhman. Stauffer. fc Mamaield l|2005j) . while the agree- 
ment is slightly worse using our new (smaller) luminos- 
ity estimate. However, the predicted luminosities are still 
well within the one a uncertainties of our measurement. 
Similarly, each of the near-IR magnitudes that we mea- 
sured for AB Dor C is consistent with the model values 
within the (large) photometric uncertainties. It is likely 
that deficiencies do exist at smaller levels in the pre- 
dicted parameters, particularly the IR m agnitudes (e.g., 
IChabrier et 31112000: iLeggett et al]l200l . 

To compare our spectral classification of AB Dor C to 
the model predictions, we must adopt a conversion of 
spectral types t o tempera t ures. Using the dwarf temper- 
ature scale_&ranjLuhma^ l)1999|) . which was a fit to data 
from lLeggett et all l|199fij) . the spectral type of M6±l for 
AB Dor C corresponds to T^e = 2840^1^^ K. This esti- 
mate does not include the uncertainty in the M dwarf 
temperature scale, which is non- negligible. An addi- 
tional uncertainty is present for AB Dor C because a 
dwarf scale, even if perfectly deter mined, may not ap- 
ply to young objects ijLuhmaiill 19991) . Nevertheless, even 
with underestimated errors, it is useful to compare our 
temperature estimate for AB Dor C to the model predic- 
tions. As shown in FigureEl our spectral classification of 
AB Dor C is con sistent with the t empe ratures predicted 
by the models of 'Chabrier et all l)2000(l when the dwarf 
scale from Luhman (1999) is adopted. 

3. INITIAL MASS FUNCTIONS OF STAR FORMING 
CLUSTERS AND THE SOLAR NEIGHBORHOOD 

The results of tests of evolutionary models at ages of 
100 Myr like the one supplied by AB Dor C do not nec- 
essarily extend to much younger ages. The uncertainties 
in the evolutionary calculations become larger at younger 
ages, and thus it would not be surprising if the models 
were robust for relatively evolved objects like AB Dor C 
that are approaching the main sequence but had signifi- 
cant errors at ages of a few million years ijBarafTe et a!.! 
120021) . Model tests that apply directly to the youngest 
ages are needed. In this section, we present such a test 
of the masses estimated with the evolutionary models of 
Haraffc et al. (1998) and Chabrier et al. ( 2000) for low- 
mass stars and brown dwarfs in star-forming regions. 

The Galactic disk appears to be populated predomi- 
nantly by stars born in embe dded clusters rather than in 
isolation ijLada fc Ladal200^ references therein). There- 
fore, unless the IMFs of embedded clusters have changed 
significantly during the lifetime of the Galaxy, the cur- 
rent generation of clusters should have the same aver- 
age IMF as the solar neighborhood. Among the best- 



studied star-forming clusters, the IMFs of stars and 
brown dwarfs inferred from a given set of evolutionary 
models exhibit relatively little variation l|Luhman et alJ 
2000). In logarithmic units where the Salpeter slope is 
1.35, the IMFs in most of these clusters rise from high 
masses down to a solar mass, rise more slowly down 
to a maximum at 0.1-0.2 Mq, and then decline into 
the substellar regime, as in the Orion Nebula Cluster 
(t ^ 0.5 Myr, Hillenbrand 119971 ILuhman et. alJl200?i 
'Hille nbrand fcC arpentcr 2000"; 'Muench et al."2002') and 
IC 348 (t ^ 2 Myr Hcrb ig 1998: Luhman_eLaL, i9983, 
120031 iMuench et 311120031) . The only clear exception is 
the IMF in the Taurus star-forming region (t ~ 1 Myr) , 
which peaks near 0.8 Mq dBriceiio ct al. 2002; L uhmanI 
120041) . but quiescent, low-density regions like Taurus pro- 
duce an insignificant number of stars compared to the 
gian t molecular clouds th at contain most embedded clus- 
ters l|Lada fc Lad3l2003() . Therefore, we C3n compare the 
model-derived IMFs of star-forming clusters to the IMF 
of the solar neighborhood to check the accuracy of mass 
estimates for young low-mass stars and brown dwarfs. 

To represent the IMF of a typical e mbedded clus- 
ter, we select the IMF of IC 348 from ILuhman et alJ 
(2003), which was derived with an H-R diagram of the 
cluster, the evolutionary models of iBaraffe et alJ l)1998f) 
and Chabrier et al. ( 20^^ . and t he temperature scale 
described in ILuhman et alJ l|2003l) . That IMF consists 
of primaries plus companions at projected separations 
greater than 300 AU, and thus more closely resembles 
an IMF of primaries rather than a single star mass func- 
tion. For the IMF of stars in the solar neighb orhood, 
we adopt the measurement by iReid et aTl l)2002|) with a 
modification to exclude companions with projected sep- 
arations less than 300 AU, making it suitable for com- 
parison to IC 348. Similar measurements of the IMF of 
field star s have been presented by Kroupa (2001) and 
IChabried l|200ll l200l . Because the mass-luminosity re- 
lation is a function of age for brown dwarfs at any age, 
and the ages of individual field brown dwarfs are un- 
known, a unique, well-sampled IMF of field brown dwarfs 
cannot be constructed. Therefore, we compare only the 
IMFs of low-mass stars between the solar neighborhood 
and IC 348. The ava ilable constraints on the substel- 
lar IMF for the field l)R,eid et all 119991: iChabrieJ l200l 
lAllen et al.ll2005l ) have been compared to measurements 
in star- forming regions in previous work ( Luhman 2004) . 

The IMFs for the solar neighborhood and the IC 348 
star-forming cluster are compared in Figure These 
IMFs are rather similar; both are roughly consistent with 
a Salpeter slope above a solar mass and are slightly ris- 
ing from a solar mass to 0.1 Mq. Previous compar- 
isons of the mass functions of young clusters and the 
solar neighborhood h ave arrived at the same conclusion 
(e.g.. lChabrieJl2(i(i^ . This agreement supports the va- 
lidity of the masses we have inferred from the models of 
Baraf fe et all (|199^ and Chabrier et al. (2000) for young 
low-mass objects. Meanwhile, C05 concluded that the 
masses of young low-mass members of clusters have been 
underestimat ed by a factor of two with the evolution- 
ary models of IChabrier et all l)2000() . To investigate this 
possibility, we have modified the IMF for IC 348 by dou- 
bling the mass estimates below 0.08 Mq. Because our 
adopted temperature scale and evolution ary models pro- 
duce accurate masses at M > 0.5 Mq ijLuhman et alJ 
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120031) , the factor by which the original masses are multi- 
plied is selected to decrease linearly with log M from two 
at 0.08 Mq to unity at 0.5 Mq. As shown in Figure [7| 
the resulting IMF differs significantly from the IMF of 
the solar neighborhood. Most notably, the relative num- 
bers of stars at 0.25 and 0.1 Mq differ by an order of 
magnitude between the two IMFs. This exercise demon- 
strates that the masses derived for the low-mass members 
of IC 348 cannot be underestimated by a factor of two if 
the IMF of Galactic star-forming clusters and the solar 
neighborhood are similar. An alternative illustration of 
this result is the following. Our use of the evolutionary 
models indicates that young objects with optical spec- 
tral types of M7 and later are brown dwarfs. Objects at 
these spectral types arc relatively rare in star- forming re- 
gions (jBriceho et al . 2002; L uhman et al... 2003: Luhman| 
l2004j) . In comparison, stars at masses of 0.1-0.2 Mq are 
the most abundant stars in the field. Therefore, members 
of star-forming regions at types of M7 and later cannot 
be low-mass stars unless the IMFs of those regions and 
the solar neighborhood are radically different. 

We have applied the above test only to the spe- 
cific methods and ni odels used in estimating masses by 
iLuhman et alJ l|2003fl . The results of this test do not ap- 
ply to the synthetic near-IR magnitudes, other sets of 
evolutionary models, or other temperature scales. 

4. DISCUSSION 

C05 found that the evolutionary models of 
iChabrier et all (jlOO^) overestimated the J and H- 
band brightnesses of AB Dor C by one magnitude and 
concluded that "the young, cool objects hitherto thought 
to be substellar in mass are therefore about twice as 
massive, which means that the frequency of brown 
dwarfs and planetary mass objects in young stellar 
clusters has been overestimated" and that "such errors 
will require serious revision of the frequency for the 
lowest-mass objects. ..at young ages". We find that this 
conclusion is oversimplified and unwarranted. First, the 
error factor in the mass estimates for young low-mass 
objects is a strong function of the stellar parameters 
considered, the choice of evolutionary models, and other 
considerations, such as the adopted conversion between 
spectral types and temperatures. For instance, even in 
the analysis of C05, the apparent error factor varied 
from 2 for J and H to only 1.3 for K. The conclusions 
of C05 included implications for measurements of sub- 
stellar IMFs in star-forming regions, but they did not 
actually test the methods predominantly used to derive 
masses in those surveys. Instead, they considered mass 
estimates based on the synthetic near-IR magnitudes, 
whereas bolometric luminosity is the most common 
parameter used in deriving masses for low-mass IMFs. 
In fact, within the age and luminosity uncertain ties 
quoted by C05, the models of lChabrier et all ("2000) are 
formally consistent with the observations of AB Dor C, 
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and so their sweeping conclusions regarding the validity 
of the models and of the measurements of substellar 
mass functions were not justified by their own data. 
Mean while, the new estimate of the a g:e for AB Dor C 
from Luh man. Stauffer. fc Mamaield ((2P05) indicates 
reasonable agreement between its predicted luminosity 
and the values measured by C05 and in this work 
fi) l2.4|l . In addition to the synthetic near-IR magnitudes, 
C05 compared the effective temperature predicted for 
AB Dor C to the spectral type of M8±l that they 
measured from IR spectroscopy, and concluded that 
the models overestimated its temperature. However, we 
measure an earlier spectral type of M6±l for AB Dor C 
from C05's spectroscopic data, which is in agreement 
with the model predictions if a dwarf temperature scale 
is adopted. Purely on an observational basis, a spectral 
type earlier than the one reported by C05 is expected 
given that AB Dor C has a mass greater than or equal 
to that of PPL 15 A and B (M6 and M7). 

We have also presented a simple test of the general 
possibility that the masses of young low-mass stars and 
brown dwarfs, particularly in star-forming regions (r ^ 
1 Myr), have been significantly underestimated in pre- 
vious studies using the models of lChabrier et al.l (|2000) . 
This test consisted of a comparison of the model-derived 
IMFs of star-forming regions and the IMF of the so- 
lar neighborhood. The agreement between these IMFs 
supports the accuracy of the masses for young low-mass 
objects when they are derived by combining measure- 
ments of spectral types and luminosities with the models 
of Baraffe et al. (1998) and Chabrier et al. j200Q) and 
the temperature scale of iLuhman et al.l l|2003fl . These 
mass estimates cannot have very large systematic errors 
(x2) unless the IMF produced by the current generation 
of star formation in the Galactic disk is significantly dif- 
ferent from the IMF of the solar neighborhood. However, 
this test does not exclude smaller systematic errors. In- 
deed, for the pri mary in a yo ung spectroscopic binary in 
Upper Scorpius, iReiners et a l. ( 2005) recently measured 
a dynamical mass of > 0.3 A/q, which is higher than 
the value of 0.2-0.3 Mq derived with the same methods 
and models used in the above IMF measurements. A 
more definitive test of those methods and models would 
be possible with a verif ication of the spect ral type of 
that binary, as noted bv IReiners et al.l (j^OOS"). Dynami- 
cal masses at lower masses are also needed beca use the 
syste matic errors could d epend on mass (|Mohantv et alJ 
,2004 iReiners et a,1 .1120(15^ . 
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TABLE 1 
Photometry of AB Dor C 



Filter 




Close et al. 


2005 






This work 








S/N 


magnitude 


F(A)/F(C) 


S/N'' 


magnitude'' 


F(A)/F(C)'' 


magnitude*^ 


F(A)/F(C)'' 


J 
H 
Ks 


10 
16 
11 


10 76+0 095 

10 04+" °'^'^ 
q 45+0.06 


150±15 
120±7.5 
80±7.5 


3.0 
4.2 
3.9 


10.72«-ji 

^- '='-0.33 


145±49 
136±33 
110±29 


10.72l»;« 

10.18t»:- 
q 7q+0.27 
^•'^-0.36 


145±64 
136±39 
110±31 



^In this measurement of AB Dor C in the long exposures, the short exposure of AB Dor A is used as the PSF. 
^Including errors from the variability of the PSF between the short and long exposures (^ V2.1.2\ . 
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Fig. 1. — The difference s of SDI images of AB Dor A at 1.625 and 1.575 /im for the two rotator an gles separate d by 33° produced by 
reprocessing the data from Clo se et all 120051) . Artificial companions at the contrast level reported by ICIose et al.l (2005) for AB Dor C 
have been inserted in the original images at four radii from AB Dor (0'.'16, 0'.'4, 0'.'6, and O'.'S). The signature of a companion produced by 
this differential imaging technique is a difference between radially shifted PSFs. The arrows indicate the location of AB Dor C reported by 
|aose et al. (2005). Neither AB Dor C nor the artificial companions with the same radius from AB Dor A are reliably detected. North is 
up and East is left in the left image. 



10 



Luhman & Potter 





Fig. 2. — Top: Near-IR images of AB Dor A (0'.'4 X 0'.'4) produced by reprocessing the data from lClose et al.1 <2005l) . Bottom: After 
applying a spatial filter to these images, AB Dor C is more noticeable. These filtered images are displayed linearly from the background 
flux at the radius of AB Dor C to the maximum flux of this companion, as represented by the grayscale bar. East is left and north is up 
in these images. 
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Fig. 3. — The spectrum of AB Dor C pro duced by reprocessing the data of IClose et'ail 1200^1^ . By comparing this spectrum to data 
for field dwarf standards fromlCushin g et alj fjDOSI) , we measure a spectral type of M6±l for AB Dor C. The standard spectra have been 
smoothed to the same resolution as the data for AB Dor C (i? ~ 1500). All spectra have been divided by polynomial fits to their continua. 
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Fig. 4. — Mfc versus J — Ks for members of the AB Dor moving group ( circZes; [Zuckerman et al."2004") and the Pleiades open cluster 
{points, T = 100-125 Myr, Mevnet et al. 1993; StaufTer et al. 1998). The mean Mk offset between the lower envelope of the Pleiades 
sequence and the observed positions of stars in the Pleiades and the AB Dor moving group are the same to within 0.05 mag, indicating 
that these two stellar populations have similar ages, just as found with Mk versus V - Ks bv rLuhman. Stauffer. fc Mamaiek L2005i'> . 
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Fig. 5. — Comparison of our measurements of the noar-IR magnitudes and bolometric luminosity of AB Dor C to the values predicted by 
the evolutionary models of Chabricr et al. 12000) for masses bracke ting its dynam ical mass of 0.09 ± 0.005 Mq Oosc et al. 2005). Using 
an age of 75-150 Myr for the AB Dor system tLuhman. Stauffer. fc Mamaiekl2005l) . the predictions are consistent with the measurements. 
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Fig. 6. — Comparison of our measurement of the spectral type of AB Dor C (M6±l) to the temp eratures predicted by the evolutionary 
models of [Chabricr ct al. (2000) for masses bracketing its dynamical mass of 0.09 i 0.005 Mq 10ose_et_alJ |2005|) . Using the dwar f 
temperature scale compiled by iLuhmani I1999f) and an age of 75-150 Myr for the AB Dor system ILuhman. Stauffer. Hi Mamaiel3l2005l) . 
the predicted temperatures are consistent with the measured spectral type. 
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Fig. 7. — Initial mass functions of stars in the solar neighbor hood (f o|)jJ^gi^^^alJ I2nn2l) and stars (solid) and brown dwarfs (dotted) in 
IC 348, which is a typical Galactic star-forming cluster (middle; Luhman ct^^200^T^ The masses in the IMF for IC 348 were derived with 
the evolutionary models of Baraffe ct al. (1998) and Chabrier et al. (2000). If the masses near and below the hydrogen burning limit in 
that IMF were significantly underestimated, then the corrected, true IMF for IC 348 (bottom) would be significantly different from that of 
the solar neighborhood. The normalization of the IMF of the solar neighborhood is arbitrary. The dashed line represents the completeness 
limit for IC 348. In the units of this diagram, the Salpeter slope is 1.35. 



